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FINITE  ELEMENT  ANALYSIS  OF  A  STATICALLY  LOADED 
ISO  TACTICAL  SHELTER 


1.  INTRODUCTION 

The  Army  has  been  using  shelters  which  are  made  from  honeycomb  panels  and  beams 
for  many  years.  The  detailed  design  of  the  shelters  has  changed  over  the  years  since 
shelters  were  designed  for  special  tasks.  Due  to  the  quantity  of  shelters  in  the  system 
and  due  to  a  desire  of  the  Army  to  have  shelters  capable  of  meeting  the  basic  requirements 
for  cargo  containers1  the  Army  has  decided  to  standardize  its  shelters. 

There  have  been  few  complete  analyses  of  shelters  in  the  past.  None  of  the  past 
studies  appeared  to  have  investigated  the  effects  of  the  panels  as  load-carrying  members 
to  determine  suitable  panel  and  beam  combinations  or  at  least  indicate  numerically  how 
the  panels  and  beams  share  the  applied  loads. 

This  study  is  a  straightforward  structural  analysis  of  a  prototype  Army  shelter 
originally  designed  to  the  International  Standard  Organization's  (ISO)  specifications  for 
dimensions  of  freight  containers.  The  objectives  of  this  study  are  to: 

(1)  Develop  a  finite  element  structural  model  for  an  ISO  Army  shelter  which  is 
2.4  m  x  2.4  m  x  6.1  m  (8  ft  x  8  ft  x  20  ft)  in  its  shipping  configuration. 

(2)  Investigate  the  load-carrying  capability  of  the  structure  by  applying  the  ISO  loads 
and  investigating  the  load  paths  determined  by  finite  element  analyses. 

(3)  Investigate  the  load-carrying  capability  of  the  structure  both  with  and  without 
panels  to  determine  the  effects  of  neglecting  panels  in  design  or  the  effects  of  using  weak 
panels  to  save  fabrication  costs. 

Section  2  of  this  report  describes  the  shelter  and  finite  element  model  chosen  for 
the  analysis.  Section  3  contains  a  description  of  the  static  loads  and  reports  the  results 
of  the  analysis  for  these  loads.  A  comparison  of  measured  data  and  computed  data  is 
made  in  Section  4.  Buckling  is  considered  for  the  stacking  load  in  Section  5  and  Section  6 
contains  the  results  and  conslusioris  obtained  from  this  numerical  study. 

2.  COMPUTER  MODEL 

a.  Description  of  Shelter 

A  general  view  of  a  3  for-1  expandable  shelter  is  shown  in  Figure  1.  The  major 
structural  components  of  the  shelter  are  aluminum  beams  and  paper  honeycomb  core 

‘ANSI  MH  5.1  —  1971,  "Basic  Requirements  for  Cargo  Containers,"  November  18,  1971. 
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sandwich  panels  with  aluminum  facings.  An  expanded  view  with  the  major  component 
breakdown  for  structural  components  in  the  shelter  is  given  in  Figure  2.  It  is  not  obvious 
in  Figure  2  that  major  load-carrying  aluminum  beams  are  built  into  the  edges  of  the  panels 
shown,  but  this  information  is  detailed  below.  Also,  there  is  a  hinge  connection  between 
the  fixed  roof  and  the  foldout  expandable  roof  which  acts  as  a  side  wall  in  the  shipping 
configuration  and  the  foldout  roof  is  not  connected  to  the  floor  when  the  shelter  is  in 
the  shipping  configuration. 

When  the  shelter  is  tested  for  static  loads  the  major  load-carrying  elements  are 
those  shown  in  Figure  3.  It  is  assumed  that  the  foldout  end  walls,  side  walls,  and  floors 
contribute  little  to  the  static  load-carrying  ability  of  the  shelter  through  their  hinge-type 
connections  to  the  structural  components  shown  in  Figure  3.  Also,  although  the  door 
frames  are  considered  as  structural  members,  the  doors  are  not,  so  this  report  contains 
the  structural  response  without  including  any  structural  interaction  which  may  occur 
between  the  door  frames  and  the  doors. 

b.  Finite  Elements  Used 

A  Government  version  (level  15.5.0)  of  the  NASTRAN  finite  element  program 
was  used  for  this  analysis.  Thus,  the  elements  chosen  for  the  analysis  were  taken  from 
the  family  of  elements  available  in  this  version  of  NASTRAN.  Since  the  shelter  is  made 
from  beams  and  sandwich  plates  the  choice  of  elements  was  dictated  by  the  problem. 
This  was  true  because  there  is  only  one  type  of  beam  element  and  only  one  type  of 
sandwich  plate  element  available  in  this  version  of  NASTRAN.  Modeling  of  hinges  required 
the  use  of  scalar  spring  elements.  Thus,  a  beam  element,  a  sandwich  plate  element,  and 
a  scalar  spring  element  was  considered  a  sufficient  variety  for  modeling  the  static  behavior 
of  the  3-for-1  shelter.  A  short  description  of  the  beam  and  plate  elements  follows.  Detailed 
information  is  available  in  the  "NASTRAN  Theoretical  Manual."2 

(1)  Beam  Element 

The  standard  straight  beam  element  was  used.  It  allows  for  a  cubic  interpolation 
ot  the  transverse  displacement  field  associated  with  bending,  a  linear  interpolation  of  the 
axial  displacement,  and  transverse  shear  and  torsion  displacements  are  treated  with  linear 
interpolation  in  the  beam.  The  element  coordinates  are  the  three  displacements  and  three 
rotations  at  each  end  of  the  beam.  Bending  is  allowed  in  two  orthogonal  planes. 

(2)  Sandwich  Plate  Element 

Triangular  Element:  The  transverse  displacement  is  interpolated  with  an 

incomplete  bicubic  polynomial  whose  coefficients  are  the  generalized  coordinates. 

2  Richard  H.  MacNeal,  "The  NASTRAN  Theoretical  Manual  (Level  15),"  COSMIC, 
University  of  Georgia,  Georgia,  April  1972. 
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3  for  1 


MAJOR  SUBASSEMBLIES 


ITEM 

DESCRIPTION 

1 

ISO  FITTING  -  LOWER  R.  H. 

(Two) 

2 

ISO  FITTING  -  LOWER  L.  H. 

(Two) 

3 

ISO  FITTING  -  UPPER  R.  H. 

(Two) 

4 

ISO  FITTING  -  UPPER  L.  H. 

(Two) 

5 

FRAME 

6 

FLOOR  PANEL 

7 

HINGED  FLOOR  PANEL 

8 

FIXED  END  PANEL  -  PERSONNEL 

DOOR 

9 

FIXED  END  PANEL  -  CARGO  DOOR 

10 

HINGED  END  PANEL 

11 

HINGED  SIDE  PANEL 

12 

FIXED  ROOF  PANEL 

13 

HINGED  ROOF  PANEL 

14 

DOOR  -  PERSONNEL 

15 

DOOR  -  CARGO 

16 

JACK 

FIGURE  2.  MAJOR  COMPONENTS  OF  A  3  FOR  1  ARMY  SHELTER 

(CONTINUED) 
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FIGURE  3.  COMPONENTS  OF  3  for-1  SHELTER  USED  FOR  STATIC  ANALYSIS 


Computation  of  stiffness  for  bending  and  transverse  shear  is  accomplished  by  expressing 
the  element's  potential  energy  as  a  function  of  bending  moments,  curvatures,  shear  forces, 
and  shear  strains.  The  material  stress-strain  laws  and  the  equilibrium  equations  are  used 
to  obtain  the  moments  and  shear  forces  in  terms  of  the  material  properties  and  strains. 
The  incomplete  bicubic  polynomial  and  the  relation  between  the  shear  strains,  slopes, 
and  rotations  are  used  to  compute  the  strains  and  curvatures  as  a  function  of  the  generalized 
coordinates.  The  stiffness  is  then  computed  relative  to  these  generalized  coordinates  by 
integration.  Finally,  the  transformation  between  the  element  coordinates  and  the 
generalized  coordinates  is  used  to  obtain  the  stiffness  matrix  in  terms  of  the  nodal 
coordinates. 

Quadrilateral  Element:  The  stiffness  matrix  is  computed  by  locating  a  median 
plane  midway  between  the  diagonals  of  the  intended  quadrilateral.  A  new  quadrilateral 
made  by  projecting  the  four  sides  of  the  original  intended  quadrilateral  (the  original  four 
points  may  not  be  in  a  plane)  onto  the  median  plane  is  used  to  compute  the  element 
stiffness  matrix.  One  diagonal  is  used  to  divide  this  new  quadrilateral  into  two  triangles. 
The  stiffness  matrices  of  these  triangles  are  then  computed,  as  stated  above,  with  respect 
to  the  degrees  of  freedom  at  the  nodes  of  the  quadrilateral.  This  process  is  repeated 
for  the  other  diagonal.  The  triangular  stiffness  matrices  are  then  weighted  by  1/2  and 
added  to  yield  the  quadrilateral  stiffness  matrix. 

c.  Finite  Elerrv.it  Mesh 

A  skelton  view  of  the  shelter  and  the  direction  of  the  axes  of  the  global 
coordinates  is  shown  in  Figure  4.  Figures  5,  6,  and  7  are  detailed  transverse,  longitudinal, 
and  plan  sectional  views  of  the  structural  components  to  be  modeled  with  the  finite 
elements.  At  first  sight  it  is  noticed  that  the  shelter  would  require  an  extremely  large 
number  of  elements,  degrees  of  freedom,  and  special  constraints  if  some  lumping  of  the 
structural  elements  is  not  done  before  discretizing  the  structure.  Thus,  the  following 
assumptions  were  made: 

(1)  Details  of  the  floor  plate  and  supporting  beam  interaction  will  not  be  sought 
in  this  analysis.  The  fact  that  the  longitudinal  spacing  between  the  transverse  beams  under 
the  floor  is  such  that  sandwich  plate  theory  is  not  even  a  valid  method  of  analysis  for 
these  plates  justifies  not  adding  degrees  of  freedom  to  more  accurately  approximate  an 
invalid  answer.  However,  it  is  hoped  (not  proved  here)  that  choosing  a  mesh  which  should 
give  a  description  of  the  floor's  bending  behavior  as  a  large  rectangular  plate  without 
the  beams  supporting  it  will  be  sufficiently  fine  for  modeling  the  floor  plate  and  beams 
as  a  major  structural  component  of  the  3-for-1  shelter.  Thus,  a  fine  mesh  for  the  floor 
plate  (to  obtain  floor  plate  data  valid  between  the  floor  beams)  was  not  made. 

(2)  Extrusions  in  the  roof,  end  walls,  and  side  walls  that  abut  and  bolt  together 
m  the  shipping  configuration  are  assumed  to  be  one  composite  beam.  This  reduces  the 
modeling  effort  by  eliminating  the  requirement  for  special  constraints  which  properly 
connect  the  degrees  of  freedom  of  two  separate  beams  to  represent  the  structural  behavior 
of  the  bolted  beams  that  may  slide  on  each  other. 
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(3)  The  hinge  connection  of  the  roof  to  the  side  wall  is  assumed  to  be  located 
at  the  intersection  of  the  midplanes  of  the  roof  and  side  wall  sandwich  plates  and  is 
not  offset  as  shown  in  Figure  5.  This  ignores  the  roof-to-wall  interaction  which  occurs 
when  the  major  (double-beam)  roof  beam  rotates  under  various  global  loading  conditions. 

(4)  The  midplanes  of  the  roof,  side  walls,  end  vails,  and  floor  are  the  only  planes 
in  which  nodal  points  are  defined,  and  the  roof,  side  walls,  etc.,  are  assumed  to  extend 
to  the  intersection  of  these  planes.  Most  of  the  column;,  and  beams  had  their  elastic 
bending  centers  close  to  these  nodes.  The  floor  stiffene  s  were  an  exception.  Their  elastic 
centers  were  thus  appropriately  offset  below  the  midplane  of  the  floor  plate.  This 
definition  of  nodal  planes  again  reduces  the  modeling  effort  and  degrees  of  freedom 
necessary  for  determining  the  global  structural  behavior  of  the  3-for-1. 

Figures  8  through  10  are  views  of  the  mesh  chosen  for  the  analysis  and  Figure  1 1 
is  a  three-dimensional  plot  of  the  finite  element  model.  The  model  consists  of  207  beams 
(10  different  properties),  54  scalar  springs,  144  quadrilateral  sandwich  plates  (3  different 
properties),  70  triangular  sandwich  plates  (3  different  properties),  and  231  nodal  points 
(1,386  degrees  of  freedom).  Thus,  the  computer  model  has  been  kept  to  a  manageable 
size  in  this  study.  Triangular  elements  were  used  (in  addition  to  quadrilateral  elements) 
to  join  domains  with  different  mesh  diameters  and  to  obtain  more  stress  output  locations 
near  the  corners. 

Figure  8  identifies  the  grid  points  and  with  reference  to  this  figure  the  following 
description  of  the  modeling  of  the  wall-to-roof  hinges  can  be  made.  In  both  cases  two 
sets  of  grid  points  are  superimposed  and  have  all  their  degrees  of  freedom  connected  by 
scalar  elements  except  the  degree  of  freedom  associated  with  the  hinge  rotation  (which 
is  left  free).  The  use  of  scalar  elements,  instead  of  constraints,  yields  data  useful  for 
approximating  the  forces  on  the  hinges.  For  this  model  both  hinges  require  rotation 
about  the  X-axis  to  be  free.  The  stiffness  of  the  scalar  elements  had  to  be  chosen  with 
care.  A  value  of  1.75T0lJN/m  was  chosen.  This  value  did  not  lead  to  ill  conditioning 
and  yielded  relative  nodal  displacements  of  the  superimposed  grid  points  of  a  value  much 
smaller  then  the  global  displacements  of  these  points.  The  first  hinge  is  represented  by 
a  one-to-one  connection  between  the  following  two  sets  of  points  (see  Figure  8):  (51, 

59,  79,  99,  119,  139,  159,  179,  199)  and  (9,  57,  77,  97,  117,  137,  157,  177,  197). 
The  second  hinge  is  given  by  a  one-to-one  connection  between  (52,  58,  78,  98,  118, 
138,  158,  178,  198)  and  (48,  56  ,  76  ,  96,  116,  136,  156,  176,  196). 

Figure  9  identifies  the  location  of  the  plate  elements  used.  The  direction  of 
the  normal  for  all  elements  in  a  given  panel  is  also  indicated  in  Figure  9.  (Note;  It 
is  important  that  the  normal  directions  be  common  in  each  panel  so  that  external  loads 
can  be  easily  applied  and  so  that  the  output  data  can  be  easily  interpreted.)  The  properties 
referred  to  in  Figure  9  are  explained  below  in  section  d,  Computation  of  Structural 
Properties. 
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Figure  10  identifies  the  location  of  all  beam  elements  used  in  the  analysis.  This 
figure  also  indicates  the  two  orthogonal  bending  planes  for  the  elements.  The  planes 
are  indicated  at  end  A  of  the  beam  element.  (Note:  Data  is  output  for  both  ends  of 
beam  elements,  hence  the  identification  of  end  A  and  end  B  in  the  form  given  in  Figure  8 
is  important  for  interpretation  of  the  output  data.)  The  properties  referred  to  in  Figure  10 
are  explained  below  in  section  d. 

The  three-dimensional  plot  given  in  Figure  11  is  useful  for  two  re  is.  First, 
this  plot  (made  by  the  computer  from  the  input  data  deck)  is  a  check  on  the  data  deck. 
Finite  element  data  decks  tend  to  be  large  and  hard  to  check.  A  plot  of  this  type  indicates 
that  at  least  the  element  and  grid  point  data  has  been  input  in  a  form  that  represents 
the  geometry  of  the  structure.  Second,  the  output  data  can  be  compared  with  the  mesh 
size  (in  more  detailed  studies)  to  evaluate  the  approximate  solution  obtained  (this  process 
requires  mesh  refinement). 

d.  Computation  of  Structural  Properties 

One  of  the  objectives  of  this  study  is  to  determine  if  the  panels  act  as  major 
load-carrying  members  (or  can  act  as  such)  when  the  shelter  is  statically  loaded  or  if 
the  frame  actually  supports  the  load  with  little  help  from  the  panels.  A  major  mechanism 
for  load  sharing  between  the  panels  and  the  frame  is  a  load  transfer  from  the  frame  to 
the  panel  in  the  plane  of  the  panel.  Thus,  good  stiffness  properties  for  the  inplane  behavior 
of  the  panel  are  required.  Since  the  mesh  was  sized  on  the  basis  of  predicting  bending 
behavior  a  computational  study  was  necessary  to  verify  that  the  mesh  chosen  for  this 
analysis  was  rot  too  stiff  to  describe  the  membrane  behavior  of  the  panels.  The  errors 
associated  with  using  a  coarse  mesh  on  membrane  problems  was  clearly  shown  by  Felippa.3 
These  errors  are  typically  associated  both  with  discretization  errors  and  with  improperly 
chosen  mesh  layouts.4'5  The  objective  of  this  computational  study  on  the  mesh  was 
to  determine  a  scale  factor  which  could  be  used  to  adjust  the  panel  stiffness  for  the 
mesh  chosen  here  to  be  near  the  panel  stiffness  of  a  properly  refined  mesh.  It  is  assumed 
that  the  magnitude  of  the  errors  associated  with  any  inconsistencies  from  the  mesh  layout 
are  smaller  than  the  errors  associated  with  modeling  the  shelter  as  an  assembly  of  beams 
and  panels. 

3C.  A.  Felippa,  "Refined  Finite  Element  Analysis  of  Linear  and  Nonlinear  Two-Dimensional 
Structures,"  Department  of  Civil  Engineering,  University  of  California,  Berkeley  Report 
66-22,  October  1966. 

4V.  Dunder  and  S.  Ridlon,  "Practical  Applications  of  Finite  Element  Method,”  Vol  104, 
No.  ST1,  ASCE  Structural  Division,  January  1978,  pp  9-21. 

J.  E.  Walz,  R.  E.  Fulton,  and  N.  J.  Cyrus,  "Accuracy  and  Convergence  of  Finite  Element 
Approximations,"  Proceedings  of  the  Second  Conference  on  Matrix  Methods  in  Structural 
Mechanics,  NTIS,  December  1968. 


The  computational  mesh  study  considered  a  2.4  m  x  2.4  m  honeycomb  panel. 
This  panel  represents  an  end  wall.  The  objective  was  to  determine  the  mesh  size  necessary 
to  obtain  good  deflection  data  for  a  racking  load  applied  in  the  form  of  a  parabolic  shear 
load  at  the  top  of  the  panel,  as  shown  in  Figure  12,  and  to  use  this  information  to 
adjust  the  inplane  stiffness  of  the  plate  elements  to  compensate  for  the  coarse  mesh.  This 
case  was  chosen  because  it  closely  represents  the  type  of  inplane  deformation  which  the 
beam  elements  can  approximate  well  but  the  inplane  constant  strain  plate  elements  cannot 
approximate  well.  The  panel  used  was  6.35  cm  thick  with  0.102  cm  thick  aluminum 
skin  and  56  kg/m3  paper  honeycomb  core.  A  computer  program  was  written  to 
automatically  generate  a  data  deck  for  NASTRAN.  The  program  was  general  enough 
so  that  any  number  of  divisions  along  the  width  and  height  could  be  chosen.  The  program 
also  computed  the  consistent  load  vector  associated  with  the  parabolic  shear  load  for  the 
given  mesh.  A  systematic  variation  of  the  mesh  was  then  used  to  obtain  the  data  shown 
in  Figure  13  and  Table  1. 

The  actual  mesh  is  nonuniform  and  has  cut-out  areas  for  doors.  A  detailed 
analysis  of  the  problems  associated  with  scaling  the  nonuniform  mesh  and  treating  strain 
singularities  at  the  cut-out  corners  is  considered  beyond  the  scope  of  this  effort.  Further, 
since  other  modeling  problems  cause  the  analysis  to  be  approximate  (lumping  the  structure 
to  plates  and  beams)  an  approximate  scale  factor  which  is  the  average  of  that  associated 
with  nearly  uniform  meshes  will  be  used.  The  end  wall  mesh  is  about  5  horizontal  elements 
by  7  vertical  elements.  This  mesh  falls  between  the  4x6  and  6x8  clases  in  Table  1. 
The  scale  factor  chosen  was  1.2  which  is  the  average  of  the  four  values  in  Table  1  for 
the  4x6  and  6x8  cases.  Thus,  the  inplane  stiffness  of  the  plates  was  multiplied  by 
1.0/1. 2  =  0.83  to  compensate  for  the  too  stiff  coarse  mesh. 

It  is  important  to  note  that  the  scaling  of  the  membrane  properties  requires 
some  extra  work  when  the  stress  data  output  by  NASTRAN  for  the  plate  elements  is 
being  used.  The  procedure  for  modifying  the  output  stress  data  can  be  explained  as 
follows.  The  output  stress  data  at  two  points  on  the  plate  cross  section  represents  the 
superposition  of  bending  stresses  and  a  membrane  stress.  The  membrane  stress  value  is 
incorrect  due  to  the  scaled  value  of  the  membrane  thickness  (scaling  was  done  to  obtain 
better  stiffness,  displacements,  and  element  loads).  The  bending  stresses  and  the  scaled 
membrane  stress  must  then  be  separated.  The  scaled  membrane  stress  should  then  be 
corrected  by  dividing  it  by  the  scale  factor  used  above  to  improve  the  element's  membrane 
stiffness.  Adding  the  correct  membrane  stress  and  the  bending  stresses  yields  the  desired 
plate  stress  data. 

There  were  three  different  types  of  sandwich  panels  used  in  the  shelter.  The 
required  structural  properties  were  the  membrane  thickness  of  the  panel,  the  transverse 
shear  thickness  of  the  panel  associated  with  the  transverse  shear  modulus,  and  the  moment 
of  inertia  of  the  panel  skins  associated  with  bending.  These  properties  were  computed 
by  hand  and  are  given  in  Table  2. 
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UNIT  SHEAR  LOAD 
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FIGURE  12.  A  2.4m  x  2.4m  PANEL  FIXED  AT  ONE  END 
WITH  SHEAR  LOAD  AT  OTHER  END. 
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DEFLECTION 


1  / N,  N  =  NUMBER  OF  ELEMENTS  PER  SIDE 


FIGURE  13  END  WALL  MESH  CONVERGENCE 
DATA  TOP  LEFT  CORNER  ONLY 
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TABLE  1.  COMPUTED  SCALE  FACTORS  FOR 
END  WALL  IN  -  PLANE  BEHAVIOR. 


MESH  SIZE 

NO.  WIDTH  x  NO.  HEIGHT 

VERTICAL 

SCALE  FACTOR 

HORIZONTAL 

SCALE 

FACTOR 

4x4 

1.32 

1.20 

8x8 

1.10 

1.08 

12  x  12 

1.06 

1.05 

16  x  16 

1.04 

1.03 

4x6 

1.27 

1.17 

6x8 

1.14 

1.10 

CONVERGED  DISPLACEMENT 

SCALE  FACTOR  =  - 

DISPLACEMENT  AT 
LISTED  MESH  SIZE 

SCALE  FACTORS  LISTED  ABOVE  ARE  AVERAGES  FOR 
POINTS  1  THROUGH  4  SHOWN  IN  FIGURE  18. 
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TABLE  2.  PHYSICAL  PROPERTIES  FOR 
PLATE  ELEMENTS 


-^^PROPERTY  NO. 

ITEM^^^-^^^ 

1 

2 

3 

MEMBRANE  THICK¬ 
NESS  SCALED  CM 

0.169 

0.169 

0.211 

MOMENT  OF  INERTIA 

CM4 

5.37 

3.46 

9.68 

TRANSVERSE  SHEAR 
THICKNESS  CM 

6.35 

5.08 

7.62 

BENDING  AND  STRET¬ 
CHING  MODULUS  PA 

6.89107 

6.89  107 

6.89  107 

MODULUS  FOR 
TRANSVERSE 

SHEAR  PA 

1.03-105 
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The  extrusions  used  in  the  shelter  being  modeled  were  special  and  their  geometric 
properties  were  not  available.  Thus  it  was  necessary  to  compute  their  properties.  Since 
there  were  ten  combinations  of  these  extrusions  and  since  they  were  not  simple  I-beams, 
a  FORTRAN  program  was  written  to  compute  their  properties.  The  program  required 
that  the  beams  be  divided  into  rectangular  sections.  The  length,  width,  location  of  center 
of  gravity  and  angle  of  inclination  was  input  for  each  rectangular  portion  of  the  beam 
and  a  program  performed  the  trivial  computations  obtaining  the  location  of  the  neutral 
axes  and  the  beam  geometric  properties.  The  ten  beam  sections  used  and  their  properties 
are  given  in  Figures  14  through  23.  Note:  In  Figures  14  through  23  the  value,  I,  indicates 
the  moment  of  inertia  about  the  1-axis.  This  value  of  inertia  is  used  for  bending  in 
plane-1  for  the  beam  as  indicated  in  Figure  10.  The  location  of  any  beam  in  the  structure 
can  be  found  by  using  Figure  10,  and  14  through  23. 

Note:  The  details  of  the  beams  and  panels  necessary  for  the  above  computations 
were  taken  from  a  preliminary  set  of  drawings.  These  drawings  represent  the  only  drawings 
available  for  the  particular  prototype  shelter  being  studied  here.  The  beam  properties 
may  be  in  slight  error  due  to  (a)  the  approximation  of  the  beam  cross  sections  as  a  group 
of  rectangles  (ignoring  curved  fillets),  and  (b)  the  lack  of  an  available  complete  data  package 
on  the  beam  cross  sections  at  the  time  when  the  computations  for  beam  properties  were 
being  verified.  Furthermore,  beam  properties  2,  3,  and  5  had  significant  cross  products 
of  inertia  which  were  not  included  in  the  analysis. 

3.  STATIC  ANALYSIS 

The  static  load  tests  given  in  the  "Specification  for  International  (ISO)  Freight 
Containers,"  reference  1,  were  divided  into  nine  separate  cases.  A  brief  description  of 
each  of  these  loading  cases  is  given  below;  for  more  detailed  information  on  the  tests 
see  either  reference  1  or  "Experimental  Measurement  of  Strain  and  Acceleration  Levels 
in  a  Rigid  Wall  Shelter  Subjected  to  Environmental  Loadings."6  All  of  these  load  cases 
were  processed  by  the  NASTRAN  finite  element  program.  The  data  obtained  on  element 
forces  was  used  to  obtain  the  results  on  load  paths  as  indicated  below.  In  order  to  clearly 
present  the  data  in  a  form  useable  to  structural  design  engineers  the  standard  format  for 
displaying  bending  moments  in  frames  used  by  civil  engineers  is  used  in  this  report.  Thus, 
the  bending  moments  are  plotted  adjacent  to  the  frame  members  with  positive  bending 
moments  inducing  tension  in  the  outer  fibers  and  positive  axial  load  for  tension  in  the 
member.  Although  shear  and  torsional  stiffnesses  were  included  in  the  frame  members 
the  loads  associated  with  shear  and  torsion  were  not  studied  here  due  to  the  effort  necessary 
to  reduce  the  data  to  graphical  form.  In  each  case  the  frame  loads  are  studied  both 
with  roof  and  wall  panels  and  without  roof  and  wall  panels.  The  data  plots  shown  below 

6F.  Barca.  'Experimental  Measurement  of  Strain  and  Acceleration  Levels  in  a  Rigid  Wall 
Shelter  Subjected  to  Evnironmental  Loadings,"  Technical  Report  TR— 79/024-AMEL,  US 
Army  Natick  Research  and  Development  Command,  Natick,  MA,  1978. 
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FIGURE  14.  BEAM  NUMBER  1  PROPERTIES 


FIGURE  15.  BEAM  NUMBER  2  PROPERTIES. 
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are  intended  to  display  both  the  distribution  of  loads  in  the  structure  as  well  as  the  change 
in  distribution  as  a  result  of  using  (or  not  using)  panels  as  structural  members.  These 
load  path  descriptions  represent  the  type  of  basic  structural  information  necessary  not 
only  for  rational  conceptual  design  studies  but  also,  more  importantly,  for  actual  designs. 
This  information  is  used  (in  normal  course  of  design)  not  only  to  size  structural  members 
but  also  to  give  the  chief  designer  the  insight  necessary  to  spot  weak  links  during  the 
detailed  design  stages  of  the  shelter  development. 

a.  Load  No.  1  -  Stacking 

The  stacking  load  represents  a  condition  that  the  shelter  is  exposed  to  in  the 
shipping  mode.  An  offset  load  (as  required  by  ISO  specifications)  is  applied  through 
the  top  ISO  fittings  and  the  floor  is  loaded  as  shown  in  Figure  24.  The  static  equivalent 
of  the  offset  load  was  computed  with  respect  to  the  nodal  locations  at  the  corners  of 
the  roof  and  these  statically  equivalent  forces  were  applied  in  the  analysis.  The  boundary 
conditions  used  in  the  analysis  are  also  shown  in  Figure  24  along  with  the  global  coordinate 
system.  The  specific  load  and  boundary  condition  data  input  can  be  found  in  the  printout 
of  the  data  deck  given  in  the  Appendix. 

Figure  25  shows  superimposed  deformed  and  undeformed  structure  plots  for 
the  beams  only.  The  maximum  deflection  was  scaled  from  0.0392  m  to  0.381  m  to 
more  clearly  display  the  form  of  the  deformation.  The  largest  deformations  occurred 
in  the  floor.  It  is  interesting  to  note  that  in  the  true  structure  the  floor  is  attached 
to  a  foldout  floor  by  a  nonstructural  hinge.  The  foldout  floor  stands  vertically  in  the 
shipping  mode  and  presents  a  very  stiff  interface  (in  the  vertical  direction)  with  the 
permanent  floor.  Failure  of  the  nonstructural  hinge  has  been  a  problem  and  the  plot 
in  Figure  25  indicates  this  could  be  due  to  the  high  flexibility  of  the  floor.  Figure  26 
is  the  deformed  body  plot  of  the  complete  finite  element  model  with  the  maximum 
deformation  again  scaled  to  0.381  m. 

It  would  not  be  practical  or  useful  to  simply  present  tables  which  give  the  output 
values  of  bending  moments,  axial  forces,  etc.,  from  the  finite  element  analysis  for  each 
load  case.  Instead,  data  was  extracted  from  several  hundred  pages  of  output  and  processed 
to  obtain  the  plots  shown  in  Figures  27  through  30.  These  plots  show  the  bending 
moment  and  axial  force  distributions  in  some  of  the  major  load-carrying  frame  members. 
The  distribution  of  these  forces  is  shown  for  two  model  cases.  As  stated  above,  one 
case  includes  the  wall  and  roof  panels  as  load-carrying  members  and  the  other  case  does 
not  include  wall  and  roof  panels.  The  benefit  of  using  panels  as  structural  members  is 
clearly  shown  in  these  plots.  The  bending  moment  in  the  vertical  frame  members  is  reduced 
significantly  except  at  the  point  of  application  of  the  moment.  A  design  using  lighter 
frame  members  could  overcome  the  high  moment  near  the  ISO  fitting  by  including  stiffener 
plates,  as  is  indeed  done  for  other  reasons  in  most  shelters.  Figure  28  indicates  that 
the  compressional  load  in  the  frame  members  is  reduced  about  30%  when  the  wall  panels 
are  used  as  structural  members. 
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FIGURE  24.  TEST  NO.  1.  STACKING 


(T)  -  Implies  bending  moments  of  this  sign 

induce  tension  in  the  outside  of  the  shelter. 
M  bending  moments 


FOR  BEAMS  IN 
FLOOR  PANEL  AND 
END  WALL  PANEL 


=  WITH  PANELS 
=  WITHOUT  PANELS 


FIGURE  77  PERSONNEL  DOOR  END.  STACKING  LOAD 

IN  PLANE  BENDING  MOMENTS  IN  OUTER  FRAME. 
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P  axial  force  in  beam,  positive  values  for  tension 


FOR  BEAMS  IN 
FLOOR  PANEL 
AND  END  WALL 
PANEL 


WITH  PANELS 
WITHOUT  PANELS 


FIGURE  28  PERSONNEL  DOOR  END.  STACKING  LOAD 
AXIAL  LOAD  IN  OUTER  FRAME. 
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FIGURE  29  SIDE  WALL  NEAREST  XZ  PLANE  STACKING  LOAD. 

IN  PLANE  BENDING  MOMENTS  IN  OUTER  FRAME 
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(T)  -  implies  bending  moments  of  this  sign  induce  tension  in  the 
bottom  of  the  beam 
M  -  bending  moment 


DISTRIBUTION  OF  BENDING  MOMENT  IN  VERTICAL 
PLANE  FOR  CENTER  CHANNEL 


- =  WITHOUT  WALL  PANELS 

-  =  WITH  WALL  PANELS 


FIGURE  30  FLOOR  BEAMS.  STACKING  LOAD. 

BENDING  MOMENTS  IN  VERTICAL  PLANE. 
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In  an  effort  to  show  that  the  same  general  results  apply  to  the  bending  moments 
in  the  plane  of  the  side  wall,  Figure  29  was  made  and  the  same  general  results  were 
found.  However,  it  should  be  noted  that  the  bending  moments  displayed  in  the  side 
wall  plane  are  much  larger  than  those  observed  in  the  end  wall  plane. 

Finally,  since  the  floor  is  uncoupled  from  the  side  walls  in  this  structure  one 
would  expect  little  change  in  the  load  distribution  iri  the  floor  when  the  wall  panels  and 
roof  panel  were  not  in  the  model.  Figure  30  shows  that  this  is  not  strictly  true.  The 
major  I-beams  under  the  floor  have  a  bending  moment  diagram  as  indicated  in  the  bottom 
set  of  curves  in  Figure  30.  The  wall  panels  tend  to  restrict  rotation  of  the  floor  at 
the  ends.  The  bending  moment  curve  with  wall  panels  acting  as  structural  members  thus 
takes  the  more  symmetric  form  as  indicated  in  the  plot.  Without  the  wall  panels  a  less 
symmetric  bending  moment  distribution  is  observed  consistent  with  the  unsymmetric 
loading.  The  top  plot  in  Figure  30  indicates  the  bending  moment  for  a  transverse  channel 
near  the  center  of  the  floor.  The  use  of  wall  and  roof  panels  has  little  effect  here.  In 
either  case  the  floor  design  does  not  appear  to  be  strongly  influenced  by  the  structural 
action  of  the  wall  panels. 

It  was  shown  above  that  the  wall  panels  significantly  reduce  the  compressive 
load  in  the  corner  posts.  To  see  the  load  distribution  in  the  wall  the  panel  stress  data 
for  the  sidewall  nearest  the  XZ  plane  was  input  to  a  FORTRAN  graphics  program  to 
obtain  Figure  31.  The  outer  skin  data  for  vertical  stresses  was  used.  The  plot  indicated 
that  the  ends  of  the  panel  near  the  corner  posts  are  in  high  compression.  This  figure 
indicates  that  the  method  of  connection  of  the  panel  to  the  frame  is  of  critical  importance. 
The  structural  design  of  the  panel-to-frame  connection  should  be  done  using  the  in-plane 
stress  data,  the  bending  moment  data,  and  the  shear  load  data  of  the  boundary  of  the 
panel. 

b.  Load  No.  2  —  Top  Lift 

The  top  lift  loading  condition  is  designed  to  test  the  shelter's  ability  to  withstand 
loads  expected  when  the  shelter  is  subjected  to  movement  by  ISO  lifting  equipment,  cranes, 
or  helicopters.  Analyses  were  made  in  this  study  for  both  a  straight  vertical  lift  at  each 
top  corner  fitting  and  for  a  central  lift.  For  the  straight  vertical  lift  data  was  obtained 
for  the  shelter  hanging  from  the  top  fittings  and  (1)  loaded  by  the  forces  generated  in 
a  vertical  gravitational  field,  (2)  loaded  by  151,460  N  uniformly  distributed  on  the  floor, 
and  (3)  a  superposition  of  (1)  and  (2).  The  central  lift  case  was  analyzed  with  the  floor 
loaded  by  151,460  N  only.  The  boundary  conditions  used  and  the  applied  loads 
(simulating  the  lifting)  are  indicated  in  Figure  32.  The  deformed  body  plot  for  the 
superimposed  straight  lift  solution  is  shown  in  Figure  33.  The  floor  deformations  are 
so  large  due  to  the  high  flexibility  of  the  floor  that  only  the  floor  appears  to  have  been 
deformed.  The  deformed  body  plot  for  the  central  lift  is  graphically  equal  to  Figure  33. 
The  maximum  deformation  in  Figure  33  was  scaled  to  be  0.381  m  from  a  computed 
value  of  0.0659  m. 
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41,302.  N  FOR  CENTRAL  LIFT  CASE 
16,521.  N  FOR  CENTRAL  LIFT  CASE 


Note,  the  gravitational  load  data  was  collected  only  for  the  lifting  tests  described 
in  this  section.  It  was  assumed  in  the  analysis  of  the  other  tests  that  the  stress  levels 
due  to  the  shelter’s  own  weight  would  not  yield  significant  design  information.  The  lifting 
tests,  however,  represented  a  complete  change  in  the  load  field  resulting  from  the  structure's 
weight  since  the  real  shelter  was  first  supported  on  the  bottom  and  then  supported  from 
the  top.  It  was  found  that  the  maximum  stresses  generated  by  the  straight  vertical  lift 
in  a  gravitational  field  with  the  floor  unloaded  were  one  order  of  magnitude  smaller  than 
the  stresses  due  to  the  straight  vertical  lift  with  only  the  floor  loaded.  Thus,  additional 
computation  designed  at  studying  the  details  of  the  element  load  fields  due  to  gravitational 
loading  was  not  in  order  for  this  design  study. 

The  stress  data  for  both  the  straight  lift  (with  gravitational  effects)  and  for  the 
central  lift  were  surveyed.  The  data  which  was  of  significant  design  value  was  almost 
identical  in  its  numerical  form  for  both  cases.  The  specific  numerical  data  which  follows 
was  abstracted  from  the  data  for  the  straight  lift  case.  Output  values  of  stress  for  the 
floor  I-beam  indicated  a  maximum  positive  extreme  fiber  bending  stress  of  2.34  •  10s  kPa. 
At  the  same  location  there  was  a  positive  axial  stress  of  0.74  •  105kPa  and  a  shear  stress 
of  3.08  •  10s  kPa.  Since  the  material  is  aluminum  (6061-T6)  with  a  yield  strength  of 
2.41  •  10s  kPa,  this  analysis  indicates  the  floor  is  not  safely  designed  for  this  type  of 
loading  (if  the  hinged  fold-out  floors  do  not  act  as  part  of  the  structure). 

The  load  distribution  data  was  investigated  for  the  corner  post  and  for  one  of 
the  longitudinal  I-beams  under  the  floor.  It  was  found  that  the  use  of  wall  panels  as 
structural  members  has  a  great  effect  on  the  magnitude  of  the  load  distribution  in  the 
shelter.  Figure  34  indicates  the  bending  moment  distribution  in  the  XZ  plane  for  the 
corner  post  nearest  the  XZ  plane  on  the  personnel  door  end.  Except  for  the  point  at 
which  the  end  wall  and  floor  interact,  the  bending  moments  in  the  corner  post  reduce 
to  insignificant  values  over  the  length  of  the  member.  Thus,  use  of  the  wall  panels  as 
load-carrying  members  and  special  design  attention  at  the  bottom  end  of  the  corner  post 
would  lead  to  optimum  design  of  the  corner  post  for  this  load  condition.  Although  the 
bending  moment  patterns  are  identical  for  the  straight  lift  and  central  lift  cases,  the  axial 
load  patterns  are  much  different.  However,  the  stresses  resulting  from  the  axial  loads 
were  small,  and  thus  a  graphical  study  of  the  difference  in  the  axial  load  paths  was  not 
made. 


Figure  35  indicates  the  bending  moments  in  the  XZ  plane  for  one  of  the 
longitudinal  I-beams  under  the  floor.  The  shape  of  the  bending  moment  curve  was  not 
changed  when  the  wall  panels  were  considered  as  part  of  the  structure  but  the  curve 
was  shifted  in  a  favorable  direction  as  the  maximum  bending  moment  was  reduced  by 
a  factor  of  about  0.6.  The  shifted  curve  takes  the  form  expected  when  the  ends  of 
the  I-beam  are  restricted  from  rotation.  In  this  case  the  wall  panels  provide  this  rotational 
interaction  at  the  ends  of  the  I-beam.  It  should  be  noted  that  this  bending  moment 
pattern  (Figure  35)  would  be  different  if  the  floor  was  structurally  attached  to  the  fold-out 
floors. 


56 


FIGURE  34  LIFTING  FROM  TOP  BENDING  MOMENTS  IN 
XZ  PLANE  FOR  CORNER  POST 


Further  investigation  of  the  stress  data  indicated  that  the  floor  plate  was  highly 
stressed.  The  longitudinal  and  transverse  stresses  were  both  large  and  near  the  yield  value 
for  the  floor  skin.  To  see  the  distribution  of  the  floor's  top  skin  stresses  in  a  form 
that  would  help  the  designer  know  where  the  worst  conditions  are,  a  three-dimensional 
plot  of  the  Hencky-von  Mises  stresses  in  the  top  skin  of  the  floor  plate  was  made.  This 
plot  is  shown  in  Figure  36.  The  highest  stressed  area  is  in  the  central  region  of  the 
floor.  The  corners  of  the  floor  plate  are  also  highly  stressed.  Thus,  in  this  case,  the 
designer  can  concentrate  his  efforts  on  satisfying  stress  requirements  in  the  central  region 
of  the  floor,  while  giving  attention  to  the  connection  methods  used  at  the  corners  of 
the  floor. 

c.  Load  No.  3  —  Bottom  Lift 

The  bottom  lift  loading  condition  is  similar  to  the  case  of  top  lifting.  As  for 
load  No.  2,  above,  the  purpose  of  load  No.  3  is  to  test  the  shelter's  ability  to  withstand 
loads  expected  from  lifting  it  from  the  bottom  ISO  fittings.  Analyses  were  made  in  this 
study  for  both  a  straight  vertical  lift  at  each  bottom  corner  fitting  and  for  a  central  lift 
by  cables  attached  between  the  bottom  ISO  fittings  on  each  side  and  one  spreader  bar 
over  the  center  of  the  shelter.  For  the  straight  vertical  lift,  data  was  obtained  for  the 
shelter  hanging  from  the  bottom  fittings  and  (1)  loaded  by  the  forces  generated  in  a 
vertical  gravitational  field,  (2)  loaded  by  151,460  N  only.  The  boundary  conditions  used 
and  the  applied  loads  are  indicated  in  Figure  37.  The  deformed  body  plot  for  the  structural 
model  with  wall  panels  was  graphically  identical  to  Figure  33  due  to  large  floor 
deformations.  Thus,  the  deformed  body  plot  for  the  case  when  the  wall  and  roof  panels 
were  not  considered  as  part  of  the  structure  was  included  as  Figure  38.  A  close  look 
at  Figure  38  indicates  that  the  end  columns  had  large  deformations  due  to  their  interaction 
with  the  floor  which  is  tending  to  rotate  at  each  end. 

As  for  load  No.  2,  above,  the  stresses  due  to  the  gravitational  loading  were 
one  order  of  magnitude  less  than  those  generated  from  loading  the  floor.  Also,  the  load 
distribution  data  of  significant  design  value  was  almost  identical  in  its  numerical  form 
for  both  the  straight  lift  and  the  central  lift  cases.  In  particular,  the  stress  data  indicated 
that  the  corner  posts  were  highly  stressed  at  the  bottom.  A  plot  of  the  bending  momenis 
in  the  corner  posts  was  made  and  the  bending  moment  distribution  was  found  *e 

identical  to  the  distribution  shown  in  Figure  34.  This  indicates  that  tfie  floor  lo. 
governed  the  moment  distribution  in  the  outer  framework  for  this  particular  structured 

The  personnel  door  frame  behaved  the  same  way  as  the  end  column.  It  was 
highly  stressed  at  the  bottom.  A  bending  moment  plot  for  the  vertical  door  frame  member 
(at  the  personnel  door  end  and  nearest  the  XZ  plane)  is  given  in  Figure  39.  The  resulting 
pattern  was  similar  to  the  pattern  for  the  corner  post.  That  is,  when  the  wall  paneis 
were  used  as  part  of  the  structure  the  panels  carried  a  significant  portion  of  the  load. 

The  stress  data  for  the  longitudinal  I  beam  under  the  floor  indicated  that  this 
member  was  highly  stressed.  The  bending  moment  pattern  for  the  I-beam  was  similar 
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(T)  -  implies  bending  moments  of  this  sign  induce 
tension  on  the  outside  of  the  shelter. 


FIGURE  39  LIFTING  FROM  BOTTOM.  BENDING 
MOMENTS  IN  XZ  PLANE  FOR 
DOOR  FRAME  MEMBER. 


to  that  shown  in  Figure  35.  The  extrusion  in  the  edge  of  the  floor  was  also  highly 
stressed.  Since  this  member  bends  with  the  I-beam  its  bending  moment  curve  is  expected 
to  be  similar  to  the  curve  for  the  I-beam.  Figure  40  indicates  that  this  is  true  and  also 
the  magnitudes  show  that  the  floor  extrusion  does  not  add  much  stiffness  to  the  floor 
structure. 


Only  one  other  area  of  the  structure  appeared  to  have  any  interesting  data  for 
this  case.  The  bottom  wall  extrusions  in  the  side  wall  were  highly  stressed  near  the  bottom 
corners  of  the  wall.  This  is  obviously  due  to  their  interaction  with  the  corner  post  (which 
is  receiving  a  rotational  load  transfer  from  the  floor). 

d.  Load  No.  4  —  Restraint 

The  restraint  test  is  a  static  test  to  check  the  design  and  construction  for  exposure 
to  horizontal  loading  conditions  expected  in  transit  of  ISO  units.  The  various  loading 
conditions  for  this  test  are  indicated  in  Figure  41.  There  are  four  possible  loadings  and 
each  was  analyzed. 

The  deformed  body  plot  for  the  entire  structure  indicated  that,  as  in  the  above 
loadings,  the  deformations  for  the  uniformly  loaded  floor  were  largest.  Figure  42  is  a 
plot  of  the  floor  plate  both  before  and  after  deformation  for  the  case  where  the  restraint 
load  is  applied  to  only  one  side  of  the  shelter  (case  No.  3  in  Figure  41).  This  deformation 
plot  indicates  that  for  this  structure  one  rnay  expect  floor  bending  to  govern  the  stress 
distribution  with  a  small  modification  from  the  presence  of  the  tension  or  compression 
load  on  the  ISO  fitting.  Figure  43  contains  interesting  information  supporting  this 
conjecture.  The  figure  contains  a  plot  of  the  axial  load  in  the  I-beam  under  the  floor 
plate  for  the  two  structural  models  (with  and  without  wall  panels)  loaded  for  case  1 
as  shown  in  Figure  41.  General  inspection  indicates  that  the  load  applied  at  the  ISO 
fitting  induces  a  tensile  load  in  the  I-beam  under  the  floor.  Figure  43  indicates  that, 
when  the  wall  panels  act  as  structural  members  the  ends  of  the  I-beam  have  zero  axial 
load.  This  appears  to  be  contradictory.  However,  careful  consideration  of  the  structure 
and  the  model  reveals  that  the  I-beam  is  offset  below  the  floor  panel.  That  is,  the  bending 
axis  of  the  I-beam  is  below  the  floor  panel  and  floor  frame  combined  structure's  neutral 
surface.  Noting  that  the  end  of  the  floor  is  more  rigidly  constrained  when  the  wall  panels 
are  part  of  the  structure,  one  would  expect  large  restraining  bending  moments  at  the 
ends  of  the  floor  (because  of  the  floor's  interaction  with  the  end  wall).  Such  restraining 
moments,  of  course,  induce  compression  below  the  neutral  surface  of  the  panel  and  I-beam 
structure  at  the  ends  of  the  floor.  When  the  end  wall  panels  are  not  part  of  the  structure, 
the  floor  panel  and  frame  combination  is  relatively  free  at  the  ends.  Bending  in  this 
case  results  in  a  lower  bending  stress  at  the  ends  of  the  floor,  thus  a  lower  counteracting 
axial  load  from  bending. 

It  is  interesting  to  note  that  although  the  floor  is  responding  in  many  ways, 
its  gross  behavior  can  be  rationally  explained  by  analogy  to  simple  bending  theory;  yet 
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(T)  -  implies  bending  moments  of  this  sign  induce  tension 
in  the  bottom  of  the  beam 

M  -  bending  moment 


=  WITH  WALL  PANELS 
=  WITHOUT  WALL  PANELS 


FIGURE  40  LIFTING  FROM  BOTTOM  BENDING  MOMENTS  IN 

XZ  PLANE  FOR  EXTRUSION  IN  EDGE  OF  FLOOR  PANEL 


65 


SIDE  WALL  NEAREST  XZ  PLANE 
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I-BEAM  UNDER  FLOOR 


-  =  WITH  WALL  PANELS 

- =  WITHOUT  WALL  PANELS 


P  =  AXIAL  LOAD  IN  I-BEAM 


FIGURE  43  RESTRAINT.  AXIAL  LOAD  IN 
I-BEAM  UNDER  FLOOR. 
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simple  bending  theory  is  inadequate  for  determining  the  numerical  data.  It  is  also  true 
that  the  modeling  of  the  structure  and  the  analysis  data  have  enlightened  the  designer 
in  the  area  of  understanding  how  this  complex  structure  works. 

Further  study  of  the  analysis  data  for  the  structural  model  containing  the  wall 
and  roof  panels  brings  out  the  following  general  design  guidance  information.  The  I-beam 
under  the  floor  was  stressed  to  about  %  of  its  yield  value.  The  stresses  for  the  entire 
shelter  were,  in  general,  small  compared  to  the  previous  cases  of  stacking  and  lifting.  The 
stress  levels  did  reach  about  %  yield  value  at  the  bottom  of  the  corner  posts  where  the 
restraint  load  was  applied.  This  area  is  stressed  from  both  the  applied  restraint  load  and 
the  uniform  loading  of  the  floor  and  is  thus  understandably  an  area  of  concern.  To 
see  how  the  side  wall  panel  responds  to  this  loading,  its  outside  skin  stresses  were  plotted 
in  Figure  44.  This  figure  serves  to  visually  demonstrate  that  the  panel  is  stressed 
significantly  only  in  the  bottom  corners  of  the  panel  near  where  the  restraint  load  is 
applied. 


When  the  wall  and  the  roof  panels  were  not  considered  as  part  of  the  structural 
model  the  stress  level  at  the  bottom  of  the  corner  posts  reduced  to  about  1/3  the  value 
it  had  when  wall  panels  were  part  of  the  structure.  However,  the  stress  level  remained 
high  in  the  remainder  of  the  corner  post.  There  was  very  little  change  in  the  stress  data 
for  the  floor  panel  and  support  frame;  that  is,  the  floor  stress  levels  were  still  far  from 
the  yield  value. 

e.  Load  No.  5  —  Racking 

The  racking  test  is  designed  to  test  the  shelter  for  its  longitudinal  and  lateral 
strength  to  resist  towing,  pushing,  and  similar  shipping  loads.  Figure  45  contains  the 
details  of  the  four  loads  and  the  boundary  conditions  used  in  this  analysis  to  represent 
the  actual  test.  This  test  is  the  first  of  the  nine  tests  studied  in  this  report  which  did 
not  have  a  uniform  load  on  the  floor  and  thus  contains  information  about  load  distributions 
in  the  end  wall  when  the  end  wall  is  not  affected  by  the  floor  bending. 

Figure  46  is  the  deformed  body  plot  for  the  entire  structure  exposed  to  the 

first  load  indicated  in  Figure  45.  The  maximum  deflection  was  scaled  to  0.254  m,  and 
thus  the  figure  is  not  as  exaggerated  as  past  figures.  A  better  indication  of  the  troubled 

area  is  given  in  Figure  47  where  a  straight-on  view  of  the  personnel  door  end  is  given 

in  the  deformed  state.  The  actual  deformation  at  the  corner  in  the  direction  of  the 
load  was  0.018  m  when  wall  and  roof  panels  acted  as  part  of  the  structure  and  0.435  m 
when  the  panels  were  not  part  of  the  structure.  Thus,  the  wall  panels  are  extremely 
important  in  keeping  the  structure  intact  when  the  shelter  is  subjected  to  racking  loads. 

After  a  survey  of  the  load  distribution  data  was  made  it  was  determined  that 
a  plot  of  the  inplane  bending  moment  distribution  in  the  outer  frame  of  the  personnel 
door  end  wall  would  yield  interesting  information.  Figure  48  contains  this  plot.  The 
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149  K N  .  i  =  LOAD  CASE  NUMBER. 


FIGURE  47.  RACKING.  LOAD  NO.  1,  DEFORMED  BODY  PLOT 

PERSONNEL  DOOR  END.  MAX.  OEF.  SCALED  TO  0.254  m. 


(T)  -  implies  bending  moments  of  this  sign  induce 
tension  in  the  outside  of  the  shelter. 

M  -  bending  moment. 


O 


WITH  PANELS  -  TOO  SMALL  TO 
BE  SEEN  ON  ABOVE  SCALE 
WITHOUT  PANELS 


FIGURE  48  RACKING.  PERSONNEL  DOOR  END.  IN  PLANE 
BENDING  MOMENTS  IN  OUTER  FRAME. 


bending  moments  in  the  frame  were  large  when  the  panels  did  not  act  as  part  of  th 
structure.  These  bending  moments  reduced  to  a  value  too  small  to  plot  on  the  sam 
scale  when  the  wall  and  roof  panels  acted  as  load-carrying  members.  Figure  48  indicate 
that  the  corners  of  the  end  wall  may  be  the  problem  area  for  design. 

The  stress  data  for  the  structural  model  with  panels  was  surveyed  and  thi 
following  general  results  obtained.  Transverse  racking  caused  the  panel  elements  in  thi 
corners  of  the  end  wall  to  have  high  inplane  principal  stresses.  These  stresses  were  als< 
high  along  the  entire  top  of  the  end  wall.  The  corner  posts  were  found  to  be  high h 
stressed  near  the  corners  of  the  end  wall,  and  the  door  frame  members  were  stressec 
to  yield  in  some  locations.  Longitudinal  racking  caused  a  behavior  in  the  side  wall  simila 
to  that  found  in  the  end  wall  above  for  transverse  racking  except  that  the  stress  magnitude 
were  well  below  yield. 

Data  for  the  model  without  panels  and  loaded  in  transverse  racking  indicatec 
high  stresses  in  the  corners  of  the  floor  panel  at  the  loaded  end  of  the  structure.  Thi 
stresses  in  the  frame  at  the  personnel  door  end  exceeded  yield  and  some  values  weri 
at  least  twice  the  ultimate  strength  of  6061— T6  aluminum.  Longitudinal  racking  alsc 
caused  high  stresses  in  the  corners  of  the  floor  plate  at  the  loaded  side  of  the  structure 
The  calculated  floor  skin  stresses  actually  reached  yield.  The  stresses  in  the  frame  member 
near  the  load  were  again  at  least  twice  yield. 

f.  Load  No.  6  —  Lashing 

The  lashing  test  loads  are  intended  to  expose  the  shelter  to  conditions  expectec 
when  the  shelter  is  secured  on  ships.  The  boundary  conditions  and  toads  used  to  obtair 
the  analysis  data  for  lashing  are  indicated  in  Figure  49.  There  are  five  separate  lashinc 
tests  in  which  opposing  compressive  and  tensile  loads  are  sequentially  applied.  Since  the 
forces  are  opposing  each  other  it  was  expected  that  the  bending  moments  and  deflection; 
found  would  be  small  relative  to  those  found  for  stacking,  lifting,  and  other  tests  in  which 
the  forces  were  not  opposing  each  other. 

The  floor  is  not  loaded  in  this  test  so  the  deformed  body  plots  are  more 
interesting  than  in  the  past  cases  where  the  floor  was  loaded  (floor  deformations  beinc 
so  large  that  they  obscured  the  remaining  deformations).  Figures  50,  51,  and  52  are 
views  of  the  personnel  door  end  as  it  is  successively  loaded  by  F1(  F3,  and  Fs  indicatec 
in  Figure  49.  The  maximum  deformations  were  scaled  from  the  computed  values  ol 
5.93  x  1CT4m,  1.24  x  1CT3m,  and  2.83  x  l(T3m  to  0.254  m,  respectively  in  Figures  5C 
to  52.  Figure  53  is  a  deformed  body  plot  of  the  entire  structure  when  it  is  loadec 
by  F8.  Figure  54  singles  out  the  side  wall  loaded  by  F8  and  shows  a  superpositior 
of  the  undeformed  and  deformed  scaled  body  plots. 

Since  the  significant  effects  of  each  loading  were  restricted  to  the  region 
surrounding  the  particular  frame  member  loaded,  it  was  decided  to  plot  the  results  ol 
several  loading  conditions  on  one  figure.  Figure  55  shows  the  axial  loads  picked  up  b\ 


FIGURE  49  TEST  NO.  6,  LASHING. 


FIGURE  50.  LASHING.  LOAD  NO.  1,  DEFORMED  BODY  PLOT 

PERSONNEL  DOOR  END.  MAX.  DEF.  SCALED  TO  0.254  m. 
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FIGURE  51.  LASHING.  LOAD  NO.  3,  DEFORMED  BODY  PLOT  PERSONNEL  DOOR  END. 
MAX.  DEF.  SCALED  TO  0.254  m. 
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FI6URE  52.  LASHING.  LOAD  NO.  5,  DEFORMED  BODY  PLOT 

PERSONNEL  DOOR  END.  MAX.  DEE.  SCALED  EO  0.254  m 
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FIGURE  55  LASHING.  AXIAL  LOADS  IN  END  WALL  FRAME. 
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the  major  frame  members  (indicated  in  the  figure)  for  load  conditions  FJ(  F3,  and  Fs. 
Figure  56  shows  similar  results  for  loadings  F8  and  F9  These  figures  indicate  that  when 
the  panels  do  not  act  as  part  of  the  structure  most  of  the  load  is  carried  by  the  beam 
between  the  opposed  loads  and  very  little,  if  any,  is  carried  away  by  the  adjacent 
framework,  the  exception  being  for  the  cases  when  frame  members  attached  to  the  floor 
plate  are  loaded  (floor  plate  was  always  considered  part  of  the  structure).  In  any  event, 
the  wall  panels  are  seen  to  reduce  the  structural  loading  by  at  least  1/3  if  these  panels 
are  properly  attached  to  the  frame  members. 

A  general  survey  of  the  stress  levels  for  load  cases  F,(  F3,  Fs,  Fg,  and  F9 
was  made  for  both  structural  models.  The  results  indicated  that  the  stress  levels  in  the 
panel's  skins  near  the  loaded  frame  member  were  about  the  same  as  the  axial  stress  levels 
in  the  frame  member.  The  results  can  be  summarized  as  follows: 

Load  Fj .  For  the  structural  model  with  wall  panels  the  end  wall  panel's  skins 
near  the  applied  load  were  stressed  to  about  1/5  yield,  as  was  the  loaded  roof  beam 
at  the  top  of  the  end  wall.  Stresses  at  the  other  locations  in  the  structure  were  negligible. 
The  model  without  the  wall  panels  had  stress  levels  in  the  roof  beam  up  to  %  yield. 

Load  F3.  The  model  with  panels  had  stress  levels  ir.  the  panel  skins  at  the 
bottom  of  the  end  wail  of  about  1/5  yield  (much  higher  applied  load  than  F,  but  a 
strong  floor  was  attached  to  help  carry  the  load).  The  transverse  beam  under  the  floor 
and  the  panel  skins  near  the  load  application  point  were  also  about  1/5  yield.  The  model 
without  panels  had  the  stresses  about  1/3  yield  in  value,  which  represents  a  50%  increase 
over  the  model  with  panels. 

Load  Fs.  The  model  with  panels  had  stresses  in  the  end  wall  panel  skins  and 
in  the  corner  post  near  the  load  application  points  of  about  14  yield  in  magnitude.  The 
panel  skin  stresses  were  only  %  yield  midway  between  the  applied  loads.  (This  may  just 
be  a  result  of  different  distances  from  the  corner  post  to  the  stress  recovery  points  in 
the  triangular  and  quadrilateral  elements.)  The  model  without  panels  had  stresses  in  the 
corner  post  of  3/4  yield.  Stresses  in  the  other  members,  of  course,  were  negligible. 

Load  F8.  The  model  with  panels  had  very  small  stresses  everywhere.  The 
highest  stresses  in  the  panel  skins  and  roof  beam  were  1/7  yield  in  value  at  the  locations 
near  the  applied  loads.  The  model  without  wall  panels  had  the  stresses  in  the  roof  beam 
increase  to  'A  yield  value,  still  relatively  small. 

Load  F9 .  This  case  is  similar  to  the  case  of  load  F3  in  that  the  floor  acts 

as  part  of  both  structural  models.  For  the  model  with  wall  panels  the  stresses  in  the 

skins  of  the  side  wall  panel  near  the  load  application  point  were  1/5  yield  as  were  the 
stresses  in  the  I-beam  under  the  floor,  the  extrusions  in  the  side  wall  panel  and  floor 

panel,  and  the  skins  of  the  floor  panel  near  the  applied  load.  The  model  without  the 

wall  panels  had  the  stresses  in  the  same  beams  and  floor  panel  skin  locations  increase 
to  about  'A  yield. 
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•IGURE  56  LASHING.  AXIAL  LOADS  IN  SIDE  WALL 
FRAME  FARTHEST  FROM  XZ  PLANE. 
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Load  No.  7  —  End  Wall  Uniform  Load 


The  end  wall  uniform  load  test  is  designed  to  test  the  shelter  for  the  cc 
expected  when  an  internal  payload  shifts  and  pushes  against  the  end  wall  durin 
at  sea.  The  loading  and  boundary  conditions  used  for  this  analysis  are  shown  in  Fi 
The  analysis  of  this  load  case  was  done  only  for  the  structural  model  with  wall  i 
panels  since  the  objective  is  to  test  the  wall.  There  was  no  door  in  the  structur 
so  the  portion  of  the  uniform  load  associated  with  the  door  was  computed  anc 
uniformly  (in  the  form  of  concentrated  loads)  around  the  door  frame. 

The  complete  structure  deformed  body  plot  is  shown  in  Figure  58.  Tf 
indicates  that  the  outward  movement  of  the  end  wall  panel  has  caused  the  flooi 
up  while  the  side  walls  and  roof  experienced  much  less  displacement  than  tl 
Figure  59  shows  the  undeformed  and  (scaled)  deformed  end  wall  plots.  Th 
indicates  that  the  highest  stresses  would  be  expected  in  the  door  frame.  Inspe 
the  stress  data  indicated  that  the  door  frame  was  the  highest  stressed  portior 
structure  but  the  stresses  were  only  on  the  order  of  1/1  Oth  the  yield  stress  for  tl 
members.  Panel  stresses  were  similarly  low.  Thus,  bending  moment,  panel  stress, 
typical  design-related  plots  were  not  made  since  it  is  felt  that  they  would  y ii 
significant  information. 

h.  Load  No.  8  —  Side  Wall  Uniform  Load 

The  purpose  of  the  side  wall  uniform  load  test  is  the  same  as  that  for 
wall  uniform  load  test.  That  is,  the  test  simulates  conditions  expected  when  an 
payload  shifts  and  pushes  against  the  side  wall  during  transit  at  sea.  The  loa< 
boundary  conditions  for  this  test  are  given  in  Figure  60.  As  in  load  No.  7  the 
was  done  only  for  the  structural  model  with  wall  and  roof  panels  since  the  < 
is  to  test  the  wall. 

The  complete  structure  deformed  body  plot  is  shown  in  Figure  61.  It  ii 
that  the  lack  of  connection  between  the  side  wall  and  the  floor  has  caused  the 
be  quite  flexible  to  the  uniform  loading  near  the  floor  edge  (otherwise  the  c 
the  panel  would  have  the  large  deflection).  It  should  also  be  noted  that  this  r 
is  conservative  since  in  the  actual  structure  the  side  wall  is  not  in  direct  cont 
the  internal  payload.  Instead,  a  fold-up  floor  hinged  to  the  main  floor  and  fold-c 
(used  for  the  expanded  structure)  are  between  the  payload  and  the  side  wall.  Thus, 
movement  of  this  load  will  be  retarded  by  the  fold-out  walls  and  hinged  fold-i 
However,  the  design  of  the  hinge  connecting  the  fold-up  floor  to  the  main  floo 
be  made  so  that  this  hinge  is  strong  enough  to  hold  the  shear  load  applied  to  i 
sliding  internal  payload.  Figure  62  is  a  detailed  superimposed  plot  for  the  und 
and  deformed  side  wall.  This  figure  indicates  the  maximum  shear  stresses  on  tl 
would  probably  occur  near  the  central  area  of  the  floor  (since  this  is  the  regio 
the  wall  flexibility  is  the  highest  and  the  loading  is  uniform). 
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This  analysis  was  for  a  3-for-1  shelter  which  has  both  side  walls  unattached 
along  the  main  floor  but  are  actually  not  directly  exposed  to  the  internal  load.  The 
2  for  1  shelter,  however,  has  one  side  wall  which  is  attached  to  the  floor  and  is  directly 
exposed  to  the  internal  payload.  This  wall  being  supported  on  all  edges  would  have  smaller 
stresses  than  those  in  the  side  wall  of  the  3-for-1  as  analyzed  above. 

Since  the  stress  levels  in  this  analysis  were  very  low  no  design  plots  of  bending 
moments,  etc.,  were  made. 

i.  Load  No.  9  —  Roof  Concentrated  Load 

The  roof  concentrated  load  test  is  derived  from  requirements  for  closed  van 
containers  and  is  designed  to  test  the  shelter  for  conditions  expected  from  small  items 
being  stacked  on  the  roof.  The  applied  load  and  boundary  conditions  used  in  this  analysis 
are  shown  in  Figure  63.  The  load  shown  was  applied  to  a  larger  area  than  indicated 
in  the  figure  so  that  a  pressure  load  could  be  applied  to  four  elements  in  the  calculation. 
This,  of  course,  would  cause  lower  stresses  to  be  computed.  Thus,  a  second  loading  model 
which  applied  the  2936  N  as  a  point  load  was  also  run  to  obtain  an  upper  bound  for 
the  computed  stresses. 

The  complete  deformed  body  plot  is  given  in  Figure  64.  Figure  65  is  a 
superimposed  plot  of  .he  undeformed  and  deformed  roof.  The  deformed  body  plots  are, 
of  course,  scaled  for  visual  effects.  The  actual  maximum  deformation  in  the  center  of 
the  roof  was  only  0.7  cm  which  is  very  small.  The  element  forces  and  stresses  were 
also  small  for  both  loading  models  mentioned  above,  and  thus  this  test  required  no  further 
study. 

4.  COMPARISON  OF  MEASURED  DATA  AND  COMPUTED  DATA 

A  series  of  tests  was  performed  on  the  shelter  being  analyzed  in  this  report.  For 
a  detailed  description  of  these  tests  see  "Experimental  Measurement  of  Strain  and 
Acceleration  Levels  in  a  Rigid  Wall  Shelter  Subjected  to  Environmental  Loadings", 
reference  6. 

A  few  comments  to  make  clear  the  meaning  of  the  computed  and  measured  data 
is  worthwhile  at  this  point.  The  real  structure  is  different  from  the  model  in  several 
ways.  The  real  structure  has  many  fastening  holes,  corner  plates  on  panels,  bolts  connecting 
two  extrusions  to  make  them  act  as  one  extrusion,  fabrication  errors,  etc.  Time  was 
not  taken  to  try  and  model  these  items  since  of  most  importance  to  the  designer  is  a 
complete  description  of  the  behavior  of  all  the  panels  and  beams  acting  together  to  resist 
the  loads.  Trying  to  model  the  details  of  fastening  holes,  etc.,  in  a  complex  structure 
of  this  type  would  complicate  the  analysis  to  an  unjustifiable  level.  The  modeling  was 
thus  directed  at  obtaining  the  gross  behavior  of  the  structure  and  so  the  comparison  of 
the  measured  and  computed  data  must  be  made  with  consideration  of  the  details  of  the 
real  structure.  A  short  discussion  of  the  measured  and  computed  data  for  each  of  the 
tests  conducted  follows. 
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FIGURE  63  TEST  NO.  9,  ROOF  CONCENTRATED  LOAD 


Stacking  Test 


One  column  at  the  personnel  door  end  was  tested  in  compression.  The  shelter 
was  laid  on  its  side  with  the  bottom  ISO  fitting  against  a  fixed  support  and  the  top 
ISO  fitting,  covered  with  a  thick  steel  plate,  had  a  hydraulic  cylinder  to  apply  the  load. 
A  computer  run  was  made  to  obtain  the  data  for  this  test.  The  stress  recovery  points 
were  at  the  extreme  fiber  points  of  the  cross  section  of  the  column.  These  stress  recovery 
points  did  not  coincide  with  the  gage  locations.  Thus,  the  element  force  data  was  used 
to  obtain  a  bending  moment  profile  and  an  axial  load  profile  for  the  column.  The  gage 
was  located  near  one  neutral  axis  of  the  column,  so  the  bending  moments  were  only 
needed  in  one  plane.  Next,  the  strains  expected  along  the  surface  of  the  column  where 
the  gages  were  mounted,  were  computed  from  the  bending  moment  and  axial  load  data. 

Figure  66  is  a  plot  of  the  computed  strain  distribution  and  the  measured  strains 
from  two  gages.  It  should  first  be  noted  that  the  corner  post  had  failed  in  a  test  done 
previous  to  this  one  and  was  not  completely  secured  to  the  rest  of  the  structure.  Thus, 
higher  measured  strain  data  would  be  expected.  The  figure  shows  good  agreement  near 
the  center  of  the  column  but  poor  agreement  near  the  bottom  of  the  column.  Also, 
poor  agreement  between  computed  and  measured  strain  data  near  the  ISO  fittings  exists 
for  other  tests  described  below.  This  indicates  that  the  modeling  near  the  corners  is 
not  sufficiently  detailed  to  accurately  compute  strain  near  the  fittings  (or  corners). 

It  should  be  noted,  however,  that  the  computation  of  strain  at  the  top  of  the 
column  indicate  a  high  strain  due  to  the  "offset"  applied  load.  In  fact,  the  computed 
value  due  to  the  "offset"  load  is  close  to  the  measured  value  found  at  the  bottom  end 
of  the  column.  A  close  inspection  of  the  column  after  the  test  indicated  that  the  fixed 
support  at  the  bottom  of  the  column  did  not  meet  the  ISO  fitting  evenly.  In  fact,  it 
was  an  "offset"  support.  This  faulty  support  may  have  induced  strains  similar  to  those 
computed  for  the  top.  This  would  then  explain  why  the  computed  and  measured  data 
do  not  agree  near  the  ISO  fitting. 

Restraint 

The  restraint  test  which  requires  the  floor  to  be  uniformly  loaded  and  the  bottom 
of  the  shelter  to  be  longitudinally  loaded  in  tension  was  made.  The  shelter's  ISO  fittings 
were  resting  on  pads  allowing  the  floor  to  bend  under  the  uniform  load  which  was  applied. 
One  end  of  the  shelter  was  restrained  from  longitudinal  translation  and  the  other  end 
fastened  to  a  cable  which  supplied  the  tensile  load.  The  stress  (strain)  recovery  points 
of  interest  were  along  the  I-beam  under  the  floor.  The  gages  were  located  at  the  neutral 
axis  of  the  I-beam.  Thus,  only  the  computed  axial  stress  data  was  needed  for  comparison 
to  the  experiment.  It  should  be  noted,  however,  that  this  axial  stress  data  for  the  I-beam 
is  a  function  of  the  floor  bending.  Also,  the  load  was  resolved  into  its  static  equivalent 
for  application  in  the  finite  element  model. 

Figure  67  is  a  plot  of  the  computed  and  measured  strains  in  the  floor  I-beam 
indicated  by  its  element  numbers.  The  computed  data  near  the  ISO  fittings  is  obviously 
not  correct.  However,  the  data  agrees  well  in  the  region  of  elements  87  to  90  which 
represents  one-half  the  length  of  the  I-beam.  The  computed  stress  data  for  the  floor 
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FIGURE  66.  STACKING  TEST.  COMPUTED  AND  MEASURED  STRAINS  IN  CORNER  POST. 
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panel,  the  extrusion  in  the  side  wall,  and  the  extrusion  in  the  edge  of  the  floor  panel 
was  surveyed  to  try  and  find  a  reason  for  the  poor  accuracy  near  the  ends  of  the  I-beam. 
No  computational  mistakes  were  found.  This  suggests  that  the  model  is  incorrect  near 
the  floor  ISO  fitting. 

The  finite  element  model  and  the  actual  shelter  were  compared  in  a  further 
attempt  to  determine  the  source  of  the  error.  It  was  noticed  that  the  assumption,  which 
states  that  the  fold-out  floors  (which  are  hinged  to  the  sides  of  the  main  floor)  do  not 
act  as  structural  members  may  be  a  poor  assumption.  The  assumption  was  originally 
made  for  the  following  reason.  The  hinge  had  failed  when  the  shelter  (loaded  with  some 
cargo)  was  moved.  This  indicates  that  the  hinge  was  not  intended  to  support  the  shear 
loads  induced  by  the  bending  of  the  floor  when  the  shelter  is  being  moved.  So  it  was 
decided  to  model  the  floor  without  including  any  stiffness  from  the  hinged  fold-out  floors. 
However,  had  the  fold-out  floor  stiffness  been  included  in  the  computations  the  axial 
compression  load  induced  by  the  bending  floor  in  the  I-beam  near  its  ends  would  not 
exist.  Thus,  the  computations  would  have  led  to  a  higher  tensile  strain  near  the  ends 
of  the  I-beam.  This  axial  strain  was  computed  and  was  found  to  be  only  —55  pm/m 
which  is  too  small  to  account  for  the  error.  Thus,  neglecting  the  hinged  fold-out  floors 
was  not  the  reason  for  the  poor  agreement  showr  in  Figure  67  near  the  ends  of  the 
I-beam. 


Further  inspection  of  the  modeling  led  to  consideration  of  the  location  of  the 
gages  with  respect  to  the  point  of  load  application  and  support.  Figure  68  indicates 
the  gage  1 1  was  located  close  to  the  fastening  point  for  the  ISO  fitting  on  the  web  of 
the  I-beam  (similarly  for  gage  15).  The  finite  element  nodes  40  and  66  are  indicated 
in  Figure  68.  The  model  assumes  that  the  axial  load  in  the  elements  between  nodes  40 
and  66  is  constant  over  the  distance  between  node  40  and  node  66  and  also  constant 
over  the  cross  section  of  that  region.  This,  of  course,  is  not  a  good  approximation  near 
the  point  where  the  load  is  applied.  Thus,  it  is  felt  that  the  coarse  mesh  near  the  corner 
is  at  least  partially  responsible  for  the  poor  agreement  between  measured  and  computed 
data  at  the  ends  of  the  I-beam. 

Lashing  —  Upper  Longitudinal 

One  upper  longitudinal  side  of  the  shelter  was  restrained  at  one  end  and  a  tensile 
load  was  applied  to  the  other  end.  Figure  69  contains  the  measured  and  computed  strains 
for  this  test.  Again,  the  agreement  is  good  except  at  the  ends.  Figure  70  is  a  detailed 
view  of  the  region  near  the  ISO  fitting  where  the  measured  strain  is  much  higher  than 
the  computed  strain.  The  corner  plate  indicated  in  the  figure  was  not  in  the  model  and 
served  to  concentrate  the  loading  near  the  end  strain  gage.  Thus,  the  measured  strains 
are  higher  than  the  computed  strains  near  the  end,  since  the  model  was  not  designed 
to  approximate  the  complicated  stress  fields  near  the  ISO  fittings. 

Lashing  —  Corner  Post  Compression 

The  corner  post  that  was  tested  for  the  stacking  load  was  tested  in  the  lashing 
compression  test.  Figure  71  indicates  the  measured  and  computed  strains.  The  central 
region,  near  gage  3,  appears  to  agree  well.  However,  the  difference  between  the  computed 
and  measured  strains  increases  near  the  ISO  fitting. 


98 


NiyyjS 


©  -  MEASURE  STRAIN. 

—  = CALCULATED  STRAIN. 

XX  *  ELEMENT  UUMSSRS. 


FIGURE  67.  RESTRAINT.  COMPUTED  AND  MEASURED  STRAINS  IN 
FLOOR  I-BEAM. 
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FIGURE  68.  RESTRAINT.  DETAILED  VIEW  OF  FLOOR  FRAME  NEAR  ISO  FITTING. 
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FIGURE  B9.  LASHING.  COMPUTED  AND  MEASURED  STRAINS  IN  ROOF  DEAM. 
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FIGURE  71.  LASHING.  COMPUTED  AND  MEASURED  STRAINS  IN  CORNER  POST. 


Lashing  —  Lower  Transverse  Compression  and  Tension 

The  tests  on  the  lower  transverse  member  had  the  most  impact  on  evaluating 
the  finite  element  modeling  and  on  determining  how  changes  in  the  modeling  can  affect 
computed  output. 

The  computed  and  measured  strains  for  the  lower  transverse  compression  and 
tension  tests  are  indicated  in  Figures  72  and  73.  The  computed  results  in  these  figures 
were  obtained  after  the  cross  sectional  area  of  beam  elements  10,  11,  167,  168,  170, 
and  171  were  changed.  The  first  comparison  of  measured  and  computed  data  was  made 
from  the  model  which  considered  the  transverse  channel  of  constant  cross  sectional  area. 
This  comparison,  not  shown  in  this  report,  was  very  poor.  The  data  indicated  that  the 
floor  plate  and  personnel  door  end  wall  were  taking  a  significant  portion  of  the  loading. 
However,  the  measured  strains  indicated  that  the  channel,  near  the  center,  was  resisting 
about  72%  of  the  applied  load  in  compression  and  94%  in  tension.  Thus,  the  structural 
stiffness  was  not  the  same  in  tension  and  compression.  Also,  the  strain  distribution  was 
constant  in  the  tension  test  but  varied  in  the  compression  test,  being  higher  near  the 
doubler  plates  in  the  compression  test  than  at  the  center.  This  means  the  structure  in 
this  area  not  only  has  different  stiffness  for  tension  and  compression  but  also  distributes 
the  load  differently  for  tension  than  for  compression. 

A  study  of  the  drawings  for  the  shelter  indicated  that  doubler  plates,  mobilizer 
fittings,  and  gusset  plates  of  high  cross  sectional  area  were  fastened  to  the  web  of  the 
transverse  channel  over  more  than  one-half  its  length.  Also,  gages  31  and  33  were  near 
the  doubler  plates.  Figure  74  is  a  detailed  view  of  the  transverse  channel.  The  figure 
contains  a  plot  of  the  total  cross  sectional  area  of  the  modified  channel  and  an  indication 
of  the  finite  element  mesh  pattern  for  the  channel  elements.  It  is  not  possible  to  represent 
the  drastic  changes  in  cross  sectional  properties  indicated  with  the  coarse  mesh  model 
used  for  this  analysis.  Thus,  a  value  of  110  cm2  was  chosen  for  elements  167,  168, 
170,  and  171  since  the  area  is  at  least  that  much  over  more  than  70%  of  those  elements. 
In  addition,  the  area  of  beam  elements  10  and  11  was  reduced  by  5.7  cm2  to  account 
for  the  modification  of  these  beams  in  the  region  of  the  personnel  door. 

The  computed  results  obtained  after  modifying  the  above  beam  properties  were 
as  expected.  The  higher  axial  stiffness  of  the  modified  channel  served  to  increase  the 
computed  strains  at  the  center  of  the  channel.  The  errors,  however,  are  still  significant 
but  the  method  of  load  distribution  was  predicted. 

5,  BUCKLING  ANALYSIS  FOR  STACKING  LOAD 

The  design  of  a  structure  should  always  include  investigations  to  determine  if  the 
structure  could  become  unstable.  The  structure  could  be  simply  underdesigned  for  the 
intended  load  or  components  of  the  structure  may  not  function  properly  causing  ?i 
unexpected  failure.  The  buckling  analysis  described  below  was  made  to  determine  if  the 
complete  structure  would  buckle  under  the  stacking  load  indicated  in  Figure  24  and  to 
determine  the  effects  that  the  wall  and  roof  panels  have  on  the  stability  of  the  structure. 
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FIGURE  72.  LASHING.  COMPUTED  AND  MEASURED  STRAINS  FOR 
LOWER  TRANSVERSE  COMPRESSION  TEST. 
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FIGURE  73.  LASHING.  COMPUTED  AND  MEASURED  STRAINS  FOR 
LOWER  TRANSVERSE  TENSION  TEST. 
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Buckling  analysis  requires  that  the  structure's  response  to  large  displacements  with 
small  strains  be  modeled.  When  this  is  done  a  formulation  can  be  made  which  determines 
a  load  scale  factor  necessary  for  a  particular  'oading  to  cause  instability,  that  is  large 
structural  displacement  without  a  further  increase  in  the  loading.  When  instability  occurs 
structural  elements  deform  to  a  configuration  not  predicted  by  the  static  analysis.  The 
determination  of  the  shape  of  this  new  displacement  field  (eigenvector),  and  the  associated 
scale  factor  (eigenvalue)  requires  the  use  of  eigenvalue  analysis  methods. 

There  are  computational  problems  involved  in  eigenanalysis  on  large  systems  which 
are  much  more  difficult  to  handle  than  those  associated  with  the  linear  small  strain  static 
analysis.  For  buckling  the  problem  is  to  find  the  lowest  eigenvalue  and  its  associated 
eigenvector.  Computing  all  eigenstates  is  prohibited  due  to  the  computer  time  required. 
The  computational  method  provided  in  the  NASTRAN  program,  which  was  used  for  this 
study,  to  determine  the  lowest  eigenstate  for  each  of  the  structural  models  is  referred 
to  as  the  "inverse  power  method  with  shifts."  As  programmed  in  NASTRAN  this  method 
requires  the  user  to  define  search  regions  in  which  the  desired  eigenvalues  are  to  be  found. 
The  user  must  also  request  the  number  of  eigenvalues  desired.  To  avoid  unnecessary 
calculations,  it  is  suggested  in  NASTRAN  literature  that  a  narrow  search  range  near  0.0 
be  chosen  for  the  eigenvalue  search  by  the  inverse  power  method.  Assuming  there  are 
no  eigenvalues  in  this  range,  the  lowest  one  outside  the  range  will  be  found.  This,  in 
turn,  should  be  the  lowest  for  the  system.  This  "narrow  band  near  0.0"  method  was 
found  to  be  very  useful. 

Tabel  3  gives  the  three  lowest  eigenvalues  for  the  structural  model  which  does  not 
have  wall  and  roof  panels.  The  table  also  includes  the  lowest  eigenvalue  associated  with 
the  complete  structure.  To  assure  that  these  values  were  correct  another  set  of 
computations  were  made  in  which  the  search  range  was  expanded  so  that  the  upper  bound 
of  the  search  range  was  slightly  less  than  the  lowest  eigenvalue  previously  predicted.  The 
same  eigenvalues  were  computed  again,  and  thus  there  are  none  below  these. 

The  table  indicates  that  if  all  the  beams  and  panels  are  acting  together  to  resist  the 
loading,  then  buckling  is  not  a  problem  since  the  structure  is  stable  until  loaded  by  6.7 
times  the  design  load.  However,  if  the  panels  are  not  properly  attached  to  the  beams 
(structure  without  wall  and  roof  panels)  then  the  critical  loading  is  only  0.52  times  the 
design  load  that  is,  buckling  will  occur  at  half  of  the  specified  load.  Thus,  this  buckling 
analysis  has  shown  that  it  is  imperative  to  use  the  panels  as  load-carrying  members  if 
the  structure  is  to  be  subjected  to  the  backing  load.  In  general,  one  could  change  the 
beam  properties  (especially  the  radius  ol  gyration)  of  the  frame  members  and  determine 
other  similar  structures  which  resist  buckiine  under  the  stacking  load  and  for  which  the 
total  frame  weight  is  less. 

It  is  interesting  to  study  the  eigenvectors  (mode  shapes)  associated  with  the  eigenvalues 
in  Table  3.  To  this  end  Figure  75  is  a  plot  of  the  first  buckling  mode  of  the  structure 
which  utilizes  the  wall  and  roof  panels  when  responding  to  the  loading.  The  plot  indicates 
that  when  all  four  corners  are  "offset"  loaded  as  indicated  in  Figure  24,  the  side  wall 
farthest  from  the  personnel  door  frame  will  bow  out  of  plane  when  the  buckling  load 
is  reached.  It  should  be  noted  that  this  behavior  of  the  complex  structure  was  not  obvious 
a  priori  and  thus  could  not  have  been  estimated  by  simpler  numerical  procedures. 
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TABLE  3.  COMPUTED  EIGENVALUES 
FOR  BUCKLING  UNDER 
STACKING  LOAD. 


MODE  NO. 

STRUCT.  NO.  1 
COMPLETE 

STRUCT.  NO.  2 
NO  PANELS 

1 

6.71 

0.525 

2 

0.681 

3 

— 

0.716 
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FIGURE  75.  BUCKLING  MODE  NO.  1  FOR  COMPLETE  STRUCTURE  SUBJECTED  TO  STACKING  LOAD. 


no 


The  first  buckling  mode  for  the  structure  without  wall  and  roof  panels  subjected 
to  the  stacking  load  is  indicated  in  Figure  76.  The  major  frame  members  buckled  into 
shapes  similar  to  those  one  would  expect  by  analogy  to  shapes  predicted  by  the  simpler 
Euler  buckling  analysis  methods  used  on  single  frame  members.  Figure  77  indicates  the 
second  mode  for  this  structure.  It  is  interesting  to  note  that  in  the  second  mode  the 
frame  members  on  the  side  farthest  from  the  personnel  door  contain  the  maximum 
amplitude  in  the  eigenvector.  This  is  similar  to  the  case  for  the  first  mode  of  the  structure 
with  wall  and  roof  panels.  However,  the  first  mode  of  the  structure  without  wall  panels 
has  its  maximum  component  unexpectedly  occurring  on  the  side  wall  nearest  to  the 
personnel  door. 

As  a  final  note,  with  respect  to  the  mode  plots  indicated  in  Figures  75  to  77,  for 
the  nonanalyst,  the  eigenvectors  computed  are  numerically  unique  only  to  a  factor  of 
—1.0.  Thus,  if  (in  testing)  one  observes  a  panel  bowing  out  instead  of  in  it  is  not  necessarily 
true  that  a  different  mode  has  been  observed. 

6.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  finite  element  model  has  been  made  for  a  3-for-1  expandable  ISO  Army  shelter. 
This  model  has  been  used  to  compute  the  gross  structural  response  of  the  frame  and 
panels  when  the  shelter  is  subjected  to  the  static  loads  specified  for  Army  ISO  Tactical 
Shelters.  The  load  paths,  which  indicate  how  the  shelter  resists  the  loading,  have  been 
displayed  and  described  in  detail  for  each  load  configuration.  These  load  paths  contain 
the  type  of  technical  information  needed  by  shelter  engineers  interested  in  refining  the 
design  of  the  3-for-1  shelter.  The  model  has  also  been  evaluated  against  measured  data 
to  determine  what  data  can  be  used  directly  from  this  type  of  finite  element  analysis 
and  what  data  must  be  determined  by  more  refined  analyses. 

The  data  presented  in  this  study  indicated  that  the  wall  panels  are  essential  structural 
members.  They  were  shown  to  share  a  major  portion  of  the  "axial"  loading  with  frame 
members.  The  wall  panels  were  also  shown  to  be  useful  for  keeping  low  bending  moments 
in  the  framework  except  near  the  ISO  fittings.  The  section  on  buckling  under  the  stacking 
load  indicated  that  the  shelter  is  much  more  stable  when  the  frame  and  panels  are  designed 
to  act  together. 

Finite  element  modeling  and  analyses  of  tactical  shelters  can  be  used  to  improve 
the  design  of  these  shelters.  As  an  example  of  how  it  can  be  used,  consider  the  main 
floor  of  the  shelter.  The  data  on  the  lift  tests  indicated  that  the  floor  was  not  designed 
well  for  these  loads.  Although  there  were  high  stresses  involved,  this  analysis  indicates 
the  basic  design  itself  is  at  fault.  The  floor  was  attached  to  fold-out  floors  along  each 
longitudinal  edge  of  the  floor  by  a  nonstructural  hinge.  The  flexibility  of  the  floor  in 
the  "Z-axis"  direction  was  shown  to  be  high  in  this  analysis.  The  "Z-axis"  flexibility 
of  the  fold-out  floors  in  the  shipping  configuration  is,  of  course,  low  since  the  panels 
are  loaded  in  plane.  This  incompatibility  of  flexibilities  is  the  design  fault.  The  result 
is  a  highly  stressed  hinge  which  may  either  fail  or  eventually  operate  improperly.  A 
suggested  fix  for  this  problem  would  be  to  connect  the  floor  to  the  side  wall  in  the 
shipping  configuration  or  to  design  the  fold-out  floor  hinges  so  that  the  fold-out  floor 
can  provide  support  for  the  main  floor. 
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FIGURE  77.  BUCKLING  MODE  NO.  2  FOR  STRUCTURE  WITH 
SUBJECTED  TO  STACKING  LOAD. 


The  data  generated  in  this  study  was  also  used  to  consider  the  feasibility  of  reducing 
the  total  mass  of  a  tactical  shelter.  The  first  suggestion  involves  the  design  of  the  edge 
extrusions  in  the  panels.  In  general,  these  extrusions  are  of  closed  cross  section.  The 
data  in  this  study  indicates  that  the  edge  extrusions  were  not  highly  stressed  in  bending 
except  near  the  ISO  fittings.  This  suggests  that  open  sections  be  considered  for  the  edge 
of  the  panels  and  the  new  designs  should  be  computationally  evaluated.  If  one  assumes 
a  20%  reduction  in  extrusion  mass  then  computations  indicate  that  1/50th  of  the  shelter's 
mass  can  be  eliminated  (approximately  47  kg).  Furthermore,  if  this  problem  were  treated 
in  detail  then  data  would  be  available  for  evaluating  the  fastening  methods  used  at  the 
boundary  of  the  panels. 

The  second  suggestion  for  reducing  the  mass  of  tactical  shelters  involves  reducing 
the  design  requirements.  If  the  shelters  could  be  given  special  treatment  in  the  shipping 
mode  then  the  column  stacking  requirements  could  be  eliminated.  A  detailed  study  of 
the  element  load  data  for  the  corner  posts  was  made  eliminating  data  from  the  stacking 
test  and  the  lashing  corner  post  compression  test.  The  maximum  bending  moments  in 
each  bending  plane  and  the  maximum  axial  load  in  the  corner  posts  were  tabulated  for 
the  remainder  of  the  test  cases.  It  was  determined  that  the  corner  post  and  floor  interaction 
put  severe  requirements  on  the  bottom  of  the  corner  post.  Assuming  this  problem  is 
resolved  by  the  doubler  plates,  gusset  plates,  etc.,  that  exist  at  the  bottoi  i  of  the  corner 
post,  the  maximum  load  data  on  the  remainder  of  the  corner  posts  was  considered.  The 
lifting  tests  and  the  transverse  racking  test  were  found  to  induce  the  highest  loads  in 
the  corner  posts.  A  conservative  redesign  of  the  corner  post  indicated  the  corner  post 
mass  could  be  cut  by  45%.  The  total  mass  reduction  for  the  four  corner  posts  was  then 
estimated  as  24  kg  (one  percent  of  total  shelter  mass).  Also,  if  the  requirements  were 
changed  this  drastically,  then  a  new  design  of  the  entire  structure  would  be  required.  How¬ 
ever  the  other  reductions  would  not  be  as  extreme. 

A  summary  of  conclusions  and  recommendations  is  as  follows: 

a.  Finite  element  modeling  of  tactical  shelters  is  useful  for  providing  technical  data 
on  load  paths  in  the  shelter. 

b.  A  coarse  finite  element  model  of  the  entire  structure  provides  load  data  valid 
everywhere  except  near  supports,  fasteners,  and  other  mechanisms  which  cause  high  stress 
gradients. 


c.  Moderate  mass  savings  may  be  obtainable  if  the  edge  extrusions  in  the  panels 
are  redesigned  for  shelters  which  have  corner  posts  separate  from  the  panels. 

d.  The  method  of  hinging  a  fold-out  panel  to  the  main  cargo  floor  is  expected 
to  place  severe  requirements  on  the  hinge  design.  If  hinge  failures  occur  it  is  suggested 
that  the  cargo  floor  be  connected  to  the  outboard  panel  in  the  shipping  mode  by  a  method 
which  reduces  the  hinge  loading. 

e.  Parameter  studies  should  be  conducted  with  the  finite  element  model.  In  these 
studies  the  design  requirements  and/or  various  frame  and  panel  structural  combinations 
could  be  considered  and  varied.  Such  a  study  would  yield  technical  data  directing  feasible 
alternative  designs. 


r 


f.  The  wall  panels  have  the  ability  to  share  significant  portions  of  the  static  loadings 
with  the  basic  framework. 

This  document  reports  research  undertaken  at  the  US  Army  Natick  Research  and 
Development  Command  and  has  been  assigned  No.  Natick/TR-79/023  in  the  series  of 
reports  approved  for  publication. 
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APPENDIX 


DATA  DECKS  FOR  FINITE  ELEMENT  MODELS 

Although  the  analysis  methods  used  to  obtain  the  data  for  the  main  body  of  this 
report  were  standard,  these  methods  are  still  abstract  to  many  engineers.  Also,  a  complete 
study  of  this  type  should  document  the  finite  element  models  in  detail.  Thus,  the  inclusion 
of  both  finite  element  models,  in  the  form  of  their  data  decks,  is  provided  here  to  help 
clarify  what  data  is  required  when  preparing  a  model  and  to  complete  the  documentation 
for  this  study. 

It  should  be  mentioned  that  not  all  the  data  is  provided.  Several  runs  were  made 
to  obtain  eigenvalue  information  for  buckling  and  special  runs  with  different  boundary 
conditions,  beam  properties,  and  applied  loads  were  made  with  temporary  data  decks. 
The  US  Customary  (inch-pound)  System  was  used  for  the  data  decks  since  the  structural 
drawings  were  in  the  US  Customary  System.  Thus,  length  is  expressed  in  inches,  force 
in  pounds,  moments  in  inch-pounds,  etc. 
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DATA  DECK  FOR  COMPLETE  SHELTER  WITH 
WALL  AND  ROOF  PANELS 
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